The structural properties of the multi-layer graphene/4H-SiC(0001) system as 

determined by Surface X-ray Diffraction 
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We present a structural analysis of the multi-layer graphene-4HSiC(0001) system using Surface 
X-Ray Reflectivity. We show for the first time that graphene films grown on the C-terminated 
(OOOl) surface have a graphene-substrate bond length that is very short (f.62A). The measured 
distance rules out a weak Van der Waals interaction to the substrate and instead indicates a strong 
bond between the first graphene layer and the bulk as predicted by ab-initio calculations. The 
measurements also indicate that multi-layer graphene grows in a near turbostratic mode on this 
surface. This result may explain the lack of a broken graphene symmetry inferred from conduction 
measurements on this system [C. Berger et al., Science 312, 1191 (2006)]. 
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INTRODUCTION 



Recent experiments have demonstrated the unique 
electronic properties of graphene sheets. P, [3, H, 0, Q 
These works point to a potential route to a new nanoelec- 
tronics paradigm based on an epitaxial graphene (EG).[H 
For the purpose of this paper we define graphene as a 
single honeycomb layer of graphite regardless of stack- 
ing order. At the moment graphene is prepared either 
by mechanical exfoliation of flakes from a bulk graphite 
sample that are subsequently deposited on an insulating 
substrate I3| or by sublimating Si from either of the 
polar faces of SiC; a process that leaves a small num- 
ber of graphene layers on the SiC surface. [H In the latter 
system, transport measurements infer that the measured 
high mobilities are limited to just a few graphene lay- 
ers (perhaps only one), that must lie near the SiC sub- 
strate. While there are similarities between the magne- 
totransport properties of exfoliated graphene and SiC- 
grown multi-layer graphene films, there are significant 
differences. [3| For instance graphene layers grown on dif- 
ferent polar faces of SiC have electron mobilities that dif- 
fer by an order of magnitude. [6] Such graphene/substrate 
specific transport properties strongly suggest that the 
substrate interaction influences the electronic properties 
of the graphene sheet. While this simple assertion may 
seem obvious, the structure and influence of the interface 
remain points of heated conjecture. One can ask if either 
the exfoliated or the SiC-grown multi-layer graphene (or 
both) are really electronically the same as a ideally iso- 
lated graphene sheet. In other words, how does the inter- 
face in both systems influence their electronic properties? 
In spite of this debate there has been no direct structural 
characterization of either the graphene-substrate inter- 



face or the graphene layers themselves in either system. 
In this work we begin to address this problem by perform- 
ing a detailed investigation of the interface structure of 
multi-layer graphene grown on the 4H-SiC(0001) surface 
using Surface X-ray Diffraction (SXD). 

Prior investigations of 6H- and 4H-SiC(0001) and 
(0001) surfaces showed that multi-layer graphene films 
can be grown on these surfaces by sublimating Si from 
SiC during heating above ~ 1200°C in ultrahigh vac- 
uum (UHV).0,Bi These studies show that multi-layer 
graphene grows epitaxially on the (0001) Si-terminated 
(Si- face) surface of SiC, while multi-layer graphene grown 
on the C-terminated (0001) (C-face) surface is rotation- 
ally disordered and under some conditions form nanocaps 
instead of a smooth film.(lo| An explanation for the 
structural differences for films grown on the two different 
faces was proposed by Forbeaux et al.[l] Their conjec- 
ture is that C-face multi-layer graphene becomes poly- 
crystalline because they have a stronger substrate-film 
bond compared to Si-face graphite. The relative bond 
strengths were inferred from K-Resolved Inverse Pho- 
toemission Spectroscopy (KRIPES),[1| and High Res- 
olution Electron Energy Loss Spectroscopy (HREELS) 
measurements. [ll| However, recent work has shown that 
the C-face multi-layer graphene is not polycrystalline and 
can be grown with domain sizes much larger than those 
grown on the Si- face. @ The improved structural order of 
C-face films correlates with magnetotransport measure- 
ments that to date find an order of magnitude improve- 
ment in electron mobilities for films grown on the C-face 
compared to Si-face films. Q Also, electronic coherence 
lengths exceeding 1 /im have been measured for multi- 
layer graphene films prepared on the C-face of SiC.[3| 
The question becomes: how can a strongly bonded C- 



face film seemingly ignore any substrate registry poten- 
tials and give rise to large free rotating films? 

Besides the question of topography differences, there 
are more fundamental questions related to electron trans- 
port in these films. For instance, how does charge 
transfer from the substrate contribute to the doping of 
graphene near the surface? One of the most important 
questions, and possibly related to the charge transport 
question, is why transport measurements on multi-layer 
graphene films grown on the (0001) C-face seem to be 
confined to just a few layers? Perhaps the most impor- 
tant question is why conduction measurements suggest 
the existence of a Berry's phase in three or more graphene 
layers. [5] Similarly, transport measurement infer'5], and 
Angle Resolved Photoemission measurements (ARPES) 
confirmfl^. the existence of a Dirac Cone in the band 
structure of multi-layer graphene films. These effects 
are not expected to occur in multi-lay er g raphene sys- 
tem with bulk graphite AB. . stacking. 13, [ij] In other 
words, experiments suggests that the major conductor is 
either an isolated single graphene sheet or possible AA.. 
stacked graphene layers. 

In this paper we begin to address these questions. 
We have performed x-ray reflectivity experiments on 
the structure of multi-layer graphene grown on the 411- 
SiC(OOOl) surface. We find that the first layer of carbon 
with an areal density of graphene sits very close to the 
last bulk SiC layer. For the C-face the graphene-bulk 
spacing is found to be 1.62 ± O.OSA. This number is con- 
sistent with recent ab-initio calculations that also indi- 
cate a covalently bonded first layer that is insulating and 
has no graphitic electronic character. [l^ We also demon- 
strate that the C-face graphene films are fiat with little 
or no corrugation in contrast with mechanically exfoli- 
ated graphite, [lij Also, by analyzing the graphite inter- 
layer spacing, we can deduce that the graphene sheets 
are stacked in a way resembling turbostratic graphite. 

These results show that films grown on the C-face of 
SiC have a strongly bonded very flat "buffer" layer. Sub- 
sequent graphene layers can be rotationally disordered 
because of the weak registry forces to this buffer layer. 
Thus the strong bonding and rotational disorder observed 
can be reconciled in a simple structural model. Most 
important, the rotational disorder and turbostratic char- 
acter of the graphite suggest that the AB.. symmetry 
of the graphite is broken leaving a graphene character 
to the films that may help to explain their conduction 
properties. 



EXPERIMENTAL 

All substrates were 4II-SiC purchased from Cree, 
Inc.JJj Prior to graphitization the 3 x 4 x 0.5mm sam- 
ples were ultrasonically cleaned in acetone and ethanol. 
Some samples were hydrogen-etched prior to graphiti- 



zation while others were not. The etching removed all 
surface scratches and left a regularly stepped surface but 
the graphite quality between etched and non-etched sam- 
ples is difficult to distinguish. Due to the difficulty of 
preparing C-face 4II-SiC samples in UHV, they were pre- 
pared by heating to 1430°C in a vacuum RF-induction 
furnace (F = 3 x 10"^ Torr) for 5-8min.[i] These pa- 
rameters produced graphitic films of 4-13 graphene lay- 
ers. The thicker films grown on the C-face refiects the 
current difficulty in producing less than 6 graphene lay- 
ers in a furnace grown environment. fil Regardless of the 
film thickness, the interface and multi-layer graphene film 
parameters measured were consistent as discussed in the 
next sections. Once samples have been graphitized, they 
remain inert allowing them to be transported into the 
Ultra High Vacuum (UHV) chamber. The x-ray scatter- 
ing experiments were performed at the Advanced Pho- 
ton Source, Argonne National Laboratory, on the 6IDC- 
^CAT UHV (F < 2 X 10-i"Torr) beam line at 16.2 keV 
photon energy. 

For all samples the number of graphene layers present 
was determined by measuring the x-ray intensity as a 
function of £ along the graphite (1,1, ^)g rod.[9| The 
notation {h,k,£)G identifies a reciprocal-space point in 
reciprocal units (r.l.u.) of the standard graphite hexag- 
onal reciprocal lattice basis vectors q = {hsLQ, khQ,£cQ, 
where a*Q — b*Q — 27r/(aGV3/2) and c*q — 2tt/cg)- The 
nominal lattice constants for graphite are qq — 2.4589A, 
cg = 6.674A.[l8'| For reflectivity data we use unsub- 
scripted reciprocal-space coordinates {h,k,£) that refer 
to the standard substrate 4H-SiC hexagonal reciprocal 
lattice units that are rotated 30° from the graphite re- 
ciprocal lattice basis. The measured lattice constants 
were: asic = 3.0802 ± 0.0006A, csic = 10.081 ± .002A 
for doped samples and asic = 3.0791 ± 0.0006A, csic = 
10.081 ± .OO2A for un-doped samples. These are within 
error bars of published values. [l9j| 



RESULTS 

To obtain detailed information about both the 
graphene films and the SiC-graphene interface, we have 
measured the surface x-ray specular refiectivity from 
graphitized 4H-SiC(0001). Specular reflectivity only de- 
pends on the momentum transfer perpendicular to the 
surface. The data is collected by integrating rocking 
curves [see Fig. Ufa)] around = for different perpen- 
dicular moment transfer vectors, q± = li^lj cgic, where 
q = k/ — ki. Since the reffectivity only depends on 
g^, the data can be analyzed using a one-dimensional 
model where all lateral information is averaged over the 
0.4 X 0.4mm x-ray beam. The scattered X-ray intensity 
/(9, g_L) is then. 



/(e,gi) - A(e,gi)|F(gi)|^ 



(1) 



where F{q±) — F{£) is the total scattering amphtude 
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FIG. 1; (a) Schematic drawing of the reflectivity geometry. 
Incident wave ki strikes sample surface at an angle O+a;. The 
diffracted wave, kf , is kept fixed at 20 from ki. q is "rocked" 
through the (001) rod by rotating the sample through an angle 
zboj. (b) Model of multi- layer graphene islands grown on a 
SiC substrate with a reconstructed SiC interface layer. For 
specular refiectivity, all n-layers islands can be represented as 
one island with a fractional surface coverage parameter, p„. 

from the 4H-SiC substrate and the graphene film over 
layer. A{Q, q±) is a. term that contains all corrections due 
to the experimental geometry. These include the Lorentz 
factor, footprint correction, effective sample area, polar- 
ization factor, and compensation for internal momentum 
transfer usin g cr itical wave vectors for both the substrate 
and the film. 2^, 2l[ F{i) is a coherent sum of the con- 
tributions from the graphene film and SiC substrate, 

F{i) = e-27s.c sin- ^i/^^{Fs,cie) + — Fg(^)}. (2) 

PSiC 

psiC and pg are the areal density of a 4H-SiC(0001) and 
a graphene (0001) plane, respectively. The weighting fac- 
tor, pG/psi=3.132, properly normalizes the scattered am- 
plitude from the graphene layer per 4H-SiC(000i) 1x1 
unit cell. The exponential term in Eq. ([21) is a rough- 
ness term that assumes csic I ^ step fluctuations in the 
SiC surface (7sic is approximately the probability of find- 
ing a step after traversing one SiC surface unit cell) . [2ll| 
The step height is measured independently from mea- 
surements of the specular rod full width at half-maximum 
(FWHM) as a function of t and gives the primary step 
height to be cs»c/2.[i,[22 

To calculate Fgicif) and Fq{€) we use a model that 
has a SiC substrate reconstruction layer and allows for 
patches of the surface to be covered with different num- 
bers of graphene layers. The schematic model of the 
graphene covered SiC is shown in Fig. [IJb) . In the model 
the SiC substrate contribution is broken into two terms: 



(i) the amplitude from a bulk terminated surface and (ii) 
the amplitude from a reconstructed SiC interface layer. 
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The first term in Eq. ([3]) is the bulk 4H-SiC structure 
factor, F^^^ if] . [23| modified by the crystal truncation 
term, (1 - e-2"^)-i,[2i] . The second term in Eq. ©, 
^iic(^)' structure factor of the reconstructed SiC 

surface. The SiC(OOOl) face is known to reconstruct into 
a 2 X 2 structure near the graphitization temperature, [isj ] 
However, the details of the reconstruction, and the nature 
of the SiC-graphene interface are not known. Although 
we cannot obtain lateral information about the inter- 
face structure from reflectivity data, the vertical shifts 
of atoms and layer density changes associated with them 
can be derived. To begin to understand this interface, 
we allow for a reconstruction by placing two SiC bilay- 
ers between the bulk and the multi-layer graphene film 
[Fig. (H^b)]. We then write the interface structure factor 
as: 



(4) 



where pj is the relative atom density for the j^^ inter- 
face layer {pj = 1 for a bulk layer) at a vertical position 
Zj.fj is the atomic form factor of C or Si. To allow for 
a reconstruction within each layer the term exp[— c^g^] 
is added. The c^'s are the rms vertical displacements 
of atoms due to a reconstruction in that layer. 26] The 
5*'' layer is added to explore the possibility of adatoms 
between the SiC and the graphene. 

To be completely general, we allow the scattered am- 
plitude, Fg{(^), from the graphene film in Eq. ^ to take 
into account the possibility of a spatial distribution of 
varying graphene layers. This is done by defining an oc- 
cupancy parameter p„ as the fractional surface area cov- 
ered by all graphene islands that are n graphene layers 
thick, pn is subject to the constraint equation = ^, 

where pq is the fraction of area that has no graphene. 
The multi-layer graphene structure factor can then be 
written in the general form: 



Faie) 




Pne 



-ql<7l/2 



E 



Do + im- l)Di m < 2 

Do + Di + {m - 2)Dg m > 2 



(5a) 



(5b) 



fc is the atomic form factor for carbon and iVmax is 
the number of graphene layers in the highest island. The 
coordinates Zm in Eq. ([5|) are the position of the to''' 
graphene atomic layer relative to the last plane of SiC 
interface atoms. Dq is the spacing between the bottom 



layer of an island and the substrate and Dq is the average 
layer spacing between graphene layers in an island (see 
Fig. H]). We have allowed the spacing between the first 
and second layer graphene, Z?i, to be different from the 
bulk to allow for differences due to a different bonding 
geometry with the substrate. Because STM studies of 
multi-layer graphene films grown on the Si-face indicate 
some buckling of the graphite layer, [i^l we also allow for a 
small vertical height distribution in each graphene layer. 
This is modelled similar to the interface relaxations in 
Eq. ([1]) by a Debye- Waller term, ac- As we will show, 
C-face graphene films show no significant buckling. 
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FIG. 2: A schematic model of multi-layer graphene grown 
on the 4H-SiC(0001) substrate. Dashed lines are the bulk 
SiC lattice planes before interface relaxation (A's). The 5"^ 
(adatom) layer is displaced Sad from the last interface layer. 
(•) are carbon atoms and (o) are silicon atoms. The shaded 
circles in the top three layers of the interface can be either 
carbon or silicon atoms depending on which of the two models, 
Carbon-corrugated or Carbon-rich, is used. 

Reflectivity data for a C-face multi-layer graphene film 
are shown in Fig. [3l The main bulk 4H-SiC peaks occur 
at £ = 4 and ^ = 8. The sharp peaks at ^ = 2, 6 and 
10 are the " quasi- forbidden" reflections of bulk SiC. [2^ 
The graphite bulk reflections are expected at £ ~ 3, 6 
and 9. While there are many variables in Eqs. ([I])-© 
that eventually must be fit, a number of the parameters 
are quite unique and insensitive to the exact structural 
model used for the SiC-graphene interface. For instance, 
because the graphite Bragg points are intense and narrow 
in i, the mean spacing between graphene layers, Dq, is 
determined with high accuracy relative to the known SiC 
lattice constant. Similarly, the graphene layer roughness 
or corrugation, ac, is determined almost solely by the 
intensity decay of the graphite Bragg points as a function 
of £. Once these nearly model-independent parameters 
are determined, they are fixed so that different structural 
models of the interface can be compared without relying 
on adjusting large numbers of parameters. 

We have tested a number of structural models for the 




FIG. 3: Specular Reflectivity vs. q± (in r.l.u.) for graphitized 
4H-SiC(0001) C-face surface with 9 graphene layers. Solid 
lines are best fits to the structural models described in the 
text. Red solid line is the best fit to the structural model 
with a smooth graphene layer (aa ~ O.OA). Blue solid line is 
the best fits with a corrugated graphite layer {ac = 0.25A). 
Magenta line is the best fit if the graphene substrate distance 
Do is reduced 10%. 



graphene/4H-SiC(0001) interface. Simple relaxations of 
the top SiC bi-layers always give poor fits to the data. 
Even attempts to make a carbon rich phase that extends 
many layers into the bulk^a model that has been pro- 
posed in the literature 0, EE ISl, was not compatible 
with the data. The best fit model is a distorted bilayer 
between the graphene and bulk SiC. A schematic of the 
model is shown in Fig. O In this model the first bilayer 
above the bulk is slightly relaxed. However, the next bi- 
layer (immediately below the graphene) has a significant 
relaxation. As we will show below, two similar versions 
of this model structure give nearly identical fits to the 
data. 

Before looking at the details of these models, we point 
out a few important model-independent parameters for 
the graphite film. First, the average graphene inter-layer 
spacmg IS found to be Dq = 3.368 ± .OOSA. This and 
other graphite film parameters are given in Table HI The 
value was determined from samples with films ranging 
from 9-13 graphene layers (averaged over the beam foot- 
print). As mentioned above the inter-layer spacing is 
nearly independent of all other fit parameters and can 
be determined with high accuracy because it is fixed by 
the £ position of the three strong graphite Bragg peaks in 
Fig. [3l The inter-layer spacing is larger than bulk crys- 
talline graphite (3.354A)[18| but smaller than the lattice 
spacmg of turbostratic graphite {Dtg = 3.440A).[33,[3lj 



TABLE I: Structural parameters for graphene grown on 4H- 
SiC(OOOl) C-Face. Parameters are defined in Fig. [2] 
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The larger spacing is due to stacking faults between ad- 
jacent layers caused by interference between tt* states 
that give rise to a larger repulsive interaction between 
adjacent graphene sheets. [30] 

Another parameter that is insensitive to the details 
of the model is, ac, in Eq. (jSap . This parameter can 
be interpreted two ways: either as a finite width of a 
graphene layer due to buckling of carbon atoms in the 
layer, or as an RMS roughness of a graphene layer due 
to vertical disorder over the coherence length of the x- 
ray beam 2/xm). We find that ac — 0.0 ± O.OSA (see 
Table [H. Because of the exponential form in Eq. ([Sa|l . a 
finite layer width manifests itself as a rapid decay in the 
graphite Bragg peak intensity at high £. This is demon- 
strated in Fig. [3] where we compare a flat graphite film to 
a film with an RMS thickness of ctg = 0.25A. The finite 
layer width severely reduces the graphite peak intensities 
at £ =6 and 9. 

Fits to the reflectivity show that two similar model 
structures for the interface region between the bulk and 
the graphene represent the experimental data equally 
well. We refer to these models as the "Carbon- 
corrugated" and "Carbon-rich" models. In both models 
the SiC bilayer immediately above the bulk in Fig. [5] re- 
mains " bulk-like" in terms of both density and inter-layer 
spacing. The two models are distinguished by the struc- 
ture of the next three layers just below the graphene film. 
Ball models of the two structures are shown in Fig. [4] 
and the detailed fitting parameters are given in Table 
im Structural values were determined for three different 
samples. The fitting parameter variations from sample to 
sample are included in the uncertainty limits of Table [TTl 

In the C-Corrugated model the last SiC bilayer is con- 
tracted inwards towards the bulk by 0.1 lA to give a 
slightly smaller Si-C bond length. In the uppermost bi- 
layer the carbon is buckled into two equal density lay- 
ers. The density of both the Si layer (pia) and the sum 
of the buckled carbon layers (pad + Pu) in this bilayer 
are each ~ 2/3 of the bulk value. It is unlikely that 
the last layer is a carbon adatom. If it were, the den- 
sity required to saturate the dangling bonds in the car- 
bon layer below would be pad — Pu/i instead of being 
equal. For this reason we refer to the model as a cor- 
rugated surface. We note that the fits are very sensitive 
to the Si density, pia, in the last bilayer. If we force the 
last bilayer to have the same Si atom density as in the 
bulk, the best fit model cannot reproduce the data. This 
is demonstrated in Fig. [5] where we show a best fit to 
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FIG. 4: Schematic ball models of bulk, C-Corrugated and C- 
Rich interface layers between the substrate and the graphene 
film. (•) are carbon atoms and (o) are silicon atoms. Hatched 
atoms are carbon atoms in the first graphene layer. Interlayer 
spacings and densities (relative to bulk SiC) are shown. 



the "C-Corrugated" model but force pia to be the bulk 
denisity. Similarly, removing the buckling in the carbon 
layer ("Smooth C-layer" model) while keeping the total 
density the same cannot reproduce the reflectivity mod- 
ulation between 0.5 < £ < 2.5 (see Fig. 

To first order the ratio of the atomic form factors for Si 
and C, fsi/ fc in Eq- Qj are determined simply by the 
ratio of their atomic numbers 14/6 — 2.33. Therefore, 
the model calculation should give a similar fit if the Si 
atoms in the top SiC bilayer are replaced by carbon atoms 
with 2.33 times the density (-- 2.33 x 0.64 = 1.49). This 
replacement gives the "C-rich" model shown in Fig. [4] 
with densities and layer spacings adjusted to give the best 
fit to the data. In Fig.|3]the best fit parameters show that 
there are two main differences between the C-Corrugated 
and C-Rich models. First, the layer spacings between 
bi-layers is considerably shorter (I.6OA) and second, the 
densities in the last layers are higher. The bilayer spacing 
measured in the C-rich model is slightly larger than the 



TABLE II: Best-fit interfacial structural parameters for graphite covered 4H-SiC(0001) (C-Face). Data for both the "C- 
Corrugated" and "C-Rich" models give nearly identical fits. All fits find Cj ~ OA for all layers. Parameters are defined in 
Fig.E] 
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bond length of diamond (1.54A). 32] The higher carbon 
layer densities have a similar significance in that they lie 
half way between the SiC density (p = 1.0) and that of 
graphene {p — 3.13). In fact, the first C-layer in the 
bilayer has a density close to the atom density of a (111) 
diamond plane, 1.51. 

While it may seem reasonable to expect that as Si sub- 
limates from the surface a carbon rich interface forms 
with some diamond- like character, we should caution 
that there are other ways to interpret these results. First 
of all, the spacing between planes in the bilayer is much 
larger, 0.63-1.03A, while in diamond they should be much 
lower, 0.5lA. The C-Rich phase is also considerably dif- 
ferent from the "extended diamond phase" proposed in 
the literature because it does not extend beyond the first 
bilayer. 0, [1^, [l^ In both models the relaxation of the 
bilayer above the bulk is small, contrary to what might 
be expected if there were significant density changes in 
that layer. These small changes from the model are 
not due to an insensitivity to either the layer spacings 
or the layer density. This is shown in Fig. [H] where we 
compare calculated best-fit reflectivities when either the 
Si density p2Si is reduced or the Si-C spacing A2Si is 
changed from the ideal value. As can be seen, inter- 
planar spacing changes of less than 5% (< 0.1 A) lead to 
obviously poor fits. Similarly, reducing the Si atom den- 
sity in this layer by more than 25% makes the fit much 
worse. Therefore, the interfacial layer does not extend 
much beyond the top-most SiC bilayer. Note also that 
the total layer density of the last three interface layers is 
p = 1.47 -I- 1.29 + 0.77 = 3.53. This density is slightly 
larger than the density of a graphene sheet (p = 3.13). 
Rather than thinking of this layer as an ideal diamond 
like layer, it may be more appropriate to view it as a 
buckled graphene sheet with a mixture of sp^ and sp^ 
bonded carbon. 

The most important finding from this work is that 
the first graphene layer sits above the last bulk carbon 
layer at a distance of Dq — 1.62 ± O.OSA. This value is, 
within error bars, insensitive to which structural model is 
used and can be determined with reasonable sensitivity as 
demonstrated in Fig. [3l The figure shows that 10% vari- 
ations in Dq from its optimal value lead to very poor fits 
to the data. The very short bond distance measured sug- 



gests that the first graphene layer is not simply bonded 
to the substrate with Van der Waal's bonds but instead 
has a much stronger interaction with the substrate. 
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FIG. 5: Specular Reflectivity for graphitized 4H-SiC(0001) 
C-face surface, (o) are the data. Red line is the best fit to the 
Carbon-corrugated top layer. The black dashed line show the 
fit for the same model if the Si layer density is fixed at the 
bulk value {pia = 1). The blue line is a fit when the carbon 
corrugation in the top layer is removed but the total density 
remains the same ("Smooth C-layer"). 



DISCUSSION 

The x-ray refiectivity data shows that the interface 
between epitaxial graphene and the 4H-SiC(0001) sub- 
strate is sharp. The interface is comprised of no more 
than 1-2 SiC bi-layers. The graphene that grows is flat 
(i.e. (Tg — OA) except for a small potential buckling of 
the first layer. There are two key structural parameters 
that deserve attention. The first is the inter-layer spac- 
ing between graphene sheets that is much larger than 
expected for AB... stacked graphene layers and points to 




FIG. 6: A comparison of the calculated reflectivity vs. q±_ (in 
r.l.u.) for different first bilayer models. Red line is the best 
fit structure with bulk bilayer parameters. Blue line is a fit 
with p2Si fixed at a value 25% less than the bulk. Magenta 
line is a best fit with the both A2Si andA2c relaxed towards 
the bulk by 5%. 



a significant stacking fault density in the film. Because 
stacking faults cause interference between tt* states in 
adjacent layers, these layers have a larger spacing. The 
mean layer spacing can, therefore, be used to estimate 
the stacking fault density, [s^l If we define the probabil- 
ity, 7, that any two adjacent sheets are faulted, then the 
inter-layer spacing will range from that of AB... stacked 
graphite (3.354A) when 7 = to that of turbostratic 
graphite (3.44A) when 7 1. In that case the aver- 
age inter-layer spacing for some finite number of stacking 
faults is approximately; [30| 

£'0 = 3.44-0.086(1-72). (6) 

Using the measured Dq — 3.368A, gives 7 = 0.4 for these 
C-face films. In other words, after every 1/(1 — 7) — 1-6 
graphene sheets, a stacking fault occurs in the film. The 
fact that there are frequent stacking faults is not sur- 
prising since there is significant rotational disorder of 
graphene layers grown on this surface. 0, Q A pair of 
graphene sheets that are rotated with respect to each 
other would lead to regions of local AB... stacking sepa- 
rated by regions with other stacking arrangements. The 
mean graphite inter-layer spacing would then be deter- 
mined by the degree of rotational disorder. Experiments 
to quantify the stacking and rotational disorder are cur- 
rently underway. 

The existence of a large stacking fault density has 
an important bearing on the results of conductivity ex- 
periments on C-face grown multi-layer graphene films. 



Magnetotransport 5], Infared Spectroscopy (IRS) [36| 
and photoemission experiments [12] indicate that multi- 
layer graphene films grown in SiC behave like isolated 
graphene sheets. In the photoemission experiments a 
clear Dirac dispersion cone is measured. The origin of 
this type of dispersion in the electronic band structure is 
the symmetry of carbon bonding in a single graphene 
sheet. AB... stacking in multi-layer graphene films 
(i.e graphite stacking) would break that symmetry and 
causes significant changes to the band structure, even 
for a few layers. 13|, |14| | In the multi-layer graphene films 
grown on the C-face of SiC, the AB... stacking disorder 
may inhibit symmetry breaking and allow sheets in the 
film to behave electronically as if they were physically 
isolated. [l3| 

In addition to the stacking fault density, the short bond 
length, Dq, between the interface and the first graphene 
sheet indicates an additional way the graphene sheets 
become isolated from the substrate. While the AB... 
stacking in bulk graphite breaks the hexagonal symme- 
try of an isolated graphene sheet, it has been assumed 
that the substrate-graphene interaction will have a sim- 
ilar effect. fl3] Indeed the short Dq bond length mea- 
sured here for the graphene/4H-SiC(0001) interface im- 
plies a strong interaction that is consistent with Photoe- 
mission results. [8] The short bond length we measure for 
the C-face films has been recently confirmed by ab-initio 
calculations. [15] Those calculations show that when a sin- 
gle graphene layer is grown on a SiC substrate the mate- 
rial remains insulating. The Dirac cone dispersion, of an 
isolated graphene sheet does not appear until a second 
graphene layer is formed. [Tsjl Therefore, the first carbon 
layer with a graphene density acts like a "buffer" layer 
between the substrate and the rest of the graphene film. 
From the structural properties of the graphene/SiC inter- 
face measured here, a model emerges that may explain 
the graphene character of these films seen in magneto- 
transport, IRS, and Photoemission measurements as well 
as in band structure calculations. In this model, one or 
two graphene layer, primarily responsible for the conduc- 
tion, lies between the "buffer" layer and the imperfectly 
stacked graphene layers above it. 

While the nature of the buffer layer is not completely 
characterized, the refiectivity data offer two possibilities. 
(1) In the C-Corrugated model the buffer layer is simply 
the first flat graphene layer above the interface. The SiC 
bilayer below this layer is relaxed with a lower density of 
atoms compared to the bulk. (2) In the C-rich structure, 
a highly buckled carbon layer, with a total carbon den- 
sity nearly equal to graphene, resides between the sub- 
strate and the rest of the film. Low Energy Electron 
Diffraction (LEED) experiments show that UHV-grown 
multi-layer graphene on the C-face surface exhibits a 2 x 2 
reconstruction. 25] Our x-ray measurements confirm that 
a 2 X 2 structure still persists even when the films are 
thick enough that LEED is no longer capable of prob- 



ing the interface. @ The long range order of the 2x2 
reconstruction is ^ 200A. This is much smaher than the 
film structural coherence length of ^ 3000A and suggests 
that the interface is never fully ordered. A possible ex- 
planation maybe that different parts of the surface are in 
different stages of graphitization. 

It is significant that the RMS layer width of the 
graphene is zero, aQ in Eq. (|5ap . ac can represent ei- 
ther a random film roughness or a RMS corrugation of 
the graphene that is commensurate with the substrate. 
Because it is zero, we can conclude that beyond the 
buffer carbon layer the graphene layers are flat and must 
be very weakly interacting with any substrate potential. 
This explains why C-face graphene films can be rota- 
tionally disordered but have large domain sizes. The en- 
ergy cost per atom to rotate a graphene sheet on a flat 
graphene substrate is very low ( < 50meV / atom) . [ssl. Is^l 
At the growth temperatures of 1400C, and given the low 
registry forces implied from these experiments, growing 
graphene sheets can rotate freely, rather than becoming 
polycrystalline as suggested by Forbeaux et al.[l] On Si- 
face multi-layer graphene films the situation is different. 
There is a (6V3 X 6V3)R300 reconstruction in the first 
2-3 graphene layers on this surface. [tI [35| The graphene 
has a nonzero corrugation of about 0.25A[23l that could 
be enough to lock the growing film into registry. Step 
boundaries or other defects in the substrate can put do- 
main boundaries in the graphene that are not easily re- 
moved by rotating large areas of the film. 

SUMMARY 

In this work we have measured a number of impor- 
tant structural parameters of multi-layer graphene grown 
on the carbon terminated 4II-SiC(000T) surface. The 
most important finding is that the separation between 
the first graphene layer and the SiC surface is very short 
(1.62 A). This distance is not much larger than a diamond 
bond length implying a very strong interaction between 
the first graphene layer and the substrate. It is consis- 
tent with recent band structure calculations that show a 
large energy gap for a single graphene layer on the 411- 
SiC(OOOl) surface.^ Subsequent graphene layers have 
an RMS corrugation (averaged over 9 layers) that is less 
than O.OSA. This suggest that the strongly bound buffer 
graphene layer shields subsequent layers from the inter- 
face corrugation potential. Therefore, unlike exfoliated 
graphene flakes deposited on Si02[lB|j the graphene lay- 
ers grown by sublimation on the C-face of SiC are very 
smooth. 

The graphene films show evidence of a large density 
of stacking faults. While the topography of these faults 
remains to be determined, it does suggest that the AB... 
stacking of bulk graphite is not present in these films. 
This may be the reason why Dirac electrons, expected 



only in isolated graphene sheets where AB... stacking 
does not break the graphene symmetry, are seen in this 
multi-graphene system. 
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